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摘 要 利用 AP-8 和 CRRESPRO 质子 辐射 带 模 式 数 据 库 , 比较 了 二 者 在 磁 赤 道 面 上 计算 结果 的 差异 并 给 出 其 
差异 原因 . 利用 AP-8 模式 数据 库 数据 , 建立 起 质子 辐射 带 函 数 形式 的 辐射 中 心 区 域 模 式 , 包括 各 能 道 全 向 微分 通 
量 峰 值 所 对 应 的 工 值 (Le) 随 能 量 E 的 变化 模式 以 及 各 能 道 全 向 微分 通 量 峰值 Jmas 随 能 量 E 的 变化 模式 . 利 
用 RBSP A 卫星 REPT 望远镜 在 磁 赤 道 面 上 的 高 能 质子 观测 数据 , 分 别 与 AP-8 模式 、CRRESPRO 模式 及 本 
文 所 得 辐射 中 心 区 域 模式 计算 结果 进行 比较 , 发 现在 78.9, 102.6 和 208 MeV 三 个 能 道上 , RBSP A 卫星 观测 所 
5 得 各 能 道 全 向 微分 通 量 明显 偏 大 , 而 Le 与 AP-8 或 本 文 辐射 中 心 模式 所 得 结果 基本 一 致 ; RBSP A 卫星 也 观测 
K. 到 CRRESPRO Quiet 模式 所 得 的 隐 性 第 二 质子 辐射 带 结构 . 
关键 词 质子 辐射 带 , 经 验 模式 , AP-8 模式 , CRRESPRO 模式 , RBSP 卫星 
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Abstract According to the databases of AP-8 and CRRESPRO proton radiation belt models, com- 
parison between both models on the magnetic equator is carried out and reasons for their differences 
are presented. The model of proton radiation belt for the peak radiation region is developed by 
using the analytic function, which includes the submodel for relationship between the L value corre- 
sponding to the peak value of omni-directional proton differential flux (Le) and the proton energy, 


and the submodel for relationship between the peak value of omni-directional proton differential flux 
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Jmax and the proton energy. Finally, the observation of energetic proton on the magnetic equator 
from RBSP A satellites REPT telescope is used to compare with the results from AP-8 model, 
CRRESPRO model and the central radiation model. It is found that the values of omni-directional 
proton differential flux for 78.9, 102.6 and 208 MeV channels from RBSP A satellite are obviously 
larger than those from the three models, but the values of Le from RBSP A satellite are very close to 
those from AP-8 model and the central radiation model, and RBSP satellite also observes the hidden 
secondary proton radiation belt presented by CRRESPRO Quiet model. 


Key words 
RBSP satellite 


0 引言 


地 球 辐射 带 又 称 Van Allen 带 , 是 指 近 地 空 间 
被 地 磁场 捕获 的 高 通 量 、 高 能 带电 粒子 集中 区 域 ， 
包括 质子 辐射 带 和 电子 辐射 带 00. 质子 辐射 带 内 
聚集 着 大 量 能 级 范围 在 0.1 到 几 百 MeV 的 高 能 质 
T. 甚至 存在 少量 达到 GeV 的 高 能 质子 D dR 
据 AP-8 模式 计算 ,大 于 10 MeV 的 高 能 质子 分 布 
在 工 为 1.14~3.7、 磁 纬度 绝对 值 约 小 于 509 的 空间 
区 域 , 在 南大 西洋 异常 区 , 质子 辐射 带 中 高 能 质子 分 
布 的 最 低 高 度 可 低 于 200 km 外， 质子 辐射 带 中 的 高 
能 质子 可 引起 航天 器 材料 的 辐射 损伤 和 电子 器 件 的 
单 粒子 效应 , 对 运行 于 其 中 的 中 低 轨 道 航天 器 在 轨 运 
行 安全 构成 威胁 , 甚至 还 威胁 到 低 轨 道 载 人 航天 器 中 
航天 员 的 生命 安全 外. 因此 , 了 解 质子 辐射 带 高 能 质 
子 的 分 布 特性 及 其 变化 规律 , 构建 质子 辐射 带 分 布 模 
式 , 对 于 减缓 或 规避 质子 辐射 带 所 造成 的 危害 起 到 至 
关 重 要 的 作用 . 

在 20 世纪 60 至 70 年 代 , 美国 NASA 依据 
1958—1970 年 大 量 卫 星 对 质子 辐射 带 的 观测 数据 ， 
构建 了 从 AP-1 到 AP-8 系列 静态 质子 辐射 带 经 
验 模式 00-7, Ep AP-8 模式 包括 太阳 活动 低 
年 MIN 和 太阳 活动 高 年 MAX MATRA, 反映 
了 太阳 辐射 对 低 高 度 质 子 辐射 带 的 有 影响， 并 且 其 适 
用 范围 较 广 , 能 量 适 用 范围 为 0.1~400 MeV, L EE 
用 范围 为 1.14~6.6, B/Bo 适用 范围 不 限 ，1994 年 
Meffert 和 Gussenhoven 利用 1990 一 1991 年 CR- 
RES 卫星 观测 数据 建立 CRRESPRO 静态 质子 辐射 
带 经 验 模式 51. 该 模式 也 包含 两 个 子 模式 , 即 地 磁 平 
静 Quiet 和 地 磁 活 路 Active 子 模式 ,反映 了 地 磁 扰动 
对 质子 辐射 带 的 影响 , 其 中 CRRESPRO Active 模式 


Proton radiation belts, Empirical model, AP-8 model, CRRESPRO model, 


还 得 到 了 第 二 质子 辐射 带 的 重要 结果 , 不 过 其 能 量 适 
用 范围 仅 为 1.1~90.4MeV, 工 值 适用 范围 为 1.0~5.5， 
B/Bo 适用 范围 为 1.0~684.6. 1998 年 Huston I P- 
fitzer 利用 1978—1995 年 TIROS-N 卫星 和 NOAA 
系列 卫星 观测 数据 , 建立 了 NOAAPRO 低 高 度 质子 
辐射 带 模式 ， 该 模式 只 适用 于 250~850km 高 度 范 
围 ， 只 能 给 出 > 16MeV, > 30MeV 和 > 80MeV 
质子 全 向 积分 通 量 , 不 过 该 模式 首次 引入 了 PFios 指 
数 因子 , 反映 太阳 活动 周 对 质子 辐射 带 的 影响 过 程 ， 
并 且 能 够 描述 Pio 与 低 高 度 三 个 能 道 质子 全 向 积 
分 通 量 之 间 的 滞后 关系 P). 1999 年 Heynderick 等 
利用 1994—1995 年 SAMPEX 卫星 PET 望远镜 
在 600km 高 度 以 下 对 18.5~500MeV 的 15 个 能 道 
高 能 质子 单 向 微分 通 量 的 观测 数据 , 建立 了 SAM- 
PEX/PET PSB97 模型 (101, 该 模式 首次 将 质子 辐射 
带 的 能 量 上 限 拓宽 到 500 MeV, 不 过 对 于 1.1~2.0 内 
不 同 的 工 值 , 实际 模式 只 能 提供 某 一 方向 上 很 罕 范 围 
内 的 投 搓 角 微分 通 量 结果 , 其 他 方向 投掷 角 微 分 通 量 
结果 只 能 通过 引入 假设 获得 . 2002 年 Huston 和 Xap- 
sos 等 利用 1978—1995 年 TIROS-N 卫星 和 NOAA 
系列 卫星 观测 数据 以 及 1990—1991 年 CRRES E 
星 观测 数据 , 构建 并 完善 了 TPM-1 (Trapped Pro- 
ton Model-1) 质子 辐射 带 经 验 模式 .该 模式 能 量 适 
用 范围 为 1.5~81.3MeV, 空间 适用 范围 为 300km 至 
地 球 静 止 同步 轨道 高 度 . 该 模式 继承 了 NOAAPRO 
模式 及 CRRESPRO 模式 的 优点 . 通过 引入 置信 度 
参数 , 使 得 该 模式 能 够 反映 未 来 太阳 活动 情况 在 不 同 
置信 水 平 下 对 质子 辐射 带 的 影响 , 其 中 不 同时 刻 所 用 
置信 和 度 值 与 对 应 太阳 活动 周 中 fio.7 峰值 通 量 以 及 
该 时 刻 在 太阳 活动 周 中 所 处 的 时 间 位 置 有 关 1719). 
自 2006 年 起 美国 多 家 机 构 联 合 利 用 以 往 多 颗 卫 星 


陶 宏 任 等 : 质子 辐射 带 辐 射 中 心 区 域 模型 


观测 数据 开发 新 一 代 的 AE-9/AP-9 辐射 带 模式 , 并 
于 2012 年 建成 AE-9/AP-9 V1.0 版 本 电子 和 质 
子 辐射 带 模式 ， 其 中 AP-9 V1.0 包括 平均 值 子 模 
Ik (Mean)、 百 分 位 子 模块 (Percentiles), 扰动 平 均 子 
模式 (Perturbed Mean) 和 蒙特 卡 罗 子 模式 (Monte 
Carlo Scenarios), 后 面 3 个 子 模式 可 设置 不 同 的 ID 
输入 , 获得 变化 结果 UU. 虽然 AP-9 V1.0 在 能 量 和 
空间 适用 范围 上 与 AP-8 模式 一 样 , 并 且 能 够 得 到 更 
复杂 的 变化 结果 , 但 该 模式 主要 是 为 了 满足 工程 设计 
需要 , 其 变化 的 输出 结果 依旧 不 能 反映 太阳 活动 周 、 
地 磁 扰 动 等 因素 的 影响 . 后 续 AE-9/AP-9 版 本 正在 
研究 中 . 

截至 目前 , 虽然 已 发 展 了 一 些 质 子 辐射 带 经 验 模 
式 , 但 对 于 不 同 能 量 在 不 同 区 域 高 能 质子 辐射 通 量 
的 计算 结果 差异 有 些 还 比较 大 ,甚至 达到 1 个 数量 
级 0- 中 . 虽然 已 了 解 到 太阳 活动 水 平 9:16717). 3 
磁 扰 动 各、 地 磁场 的 长 期 衰减 和 西 漂 167201. HERO 
气 数 密度 等 因素 对 质子 辐射 带 的 影响 DU, 但 仍 缺 乏 
能 够 描述 太阳 活动 、 地 磁 扰 动 等 因素 影响 的 动态 变 
化 的 整个 质子 辐射 带 经 验 模式 ， 其 主要 原因 是 卫星 
观测 数据 质量 (例如 数据 污染 、 数 据 定 标 等 问题 ) 及 
网 格 化 数据 库 建 模 方 法 的 局 限 性 . 

考虑 到 以 往 质子 辐射 带 观 测 数据 的 质量 问题 , 本 
文 将 利用 AP-8 模式 数据 库 , 建立 函数 形式 的 质子 辐 
射 带 辐 射 中 心 区域 模 式 , 为 构建 函数 形式 的 动态 质子 
辐射 带 经 验 模式 奠定 基础 .同时 , 将 利用 目前 RB- 
SP 卫星 REPT 望远镜 高 能 质子 观测 结果 与 AP-8 
模式 、CRRESPRO 模式 及 质子 辐射 带 辐射 中 心 区 域 
模式 计算 结果 进行 对 比 . 


1 数据 介绍 


本 文 所 用 数据 来 自 AP-8 模式 数据 库 、CRRE- 
SPRO 模式 数据 库 及 RBSP 卫星 上 相对 论 电 子 质子 
望远镜 (REPT) 高 能 质子 观测 数据 . AP-8 模式 数据 
库 数 据 用 于 辐射 带 辐射 中 心 区 域 建 模 , 其 余 两 个 用 于 
与 AP-8 模式 数据 库 结果 或 辐射 中 心 区 域 模式 结果 
的 比较 . 

1.1 AP-8 模式 数据 库 

AP-8 模式 包含 了 AP-8 MIN 太阳 活动 低 年 子 模 

式 和 AP-8 MAX 太阳 活动 高 年 子 模 式 . 这 两 个 子 模 


*http://craterre.onecert.fr/prbem /irbem/description.html 
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式 都 是 静态 平均 模式 , 均 采 用 后 台 网 格 化 数据 库 的 形 
式 建 模 , 通过 查询 数据 库 及 相应 的 插值 方法 获得 质子 
辐射 带 中 各 能 量 的 高 能 质子 辐射 情况 . 

AP-8 模式 后 台 网 格 数据 库 采 用 (E, L, B/ Bo) 
网 格 点 ,保存 各 网 格 点 上 相应 能 量 的 全 向 积分 通 量 . 
AP-8 模式 后 台 网 格 数据 库 E 网 格 点 的 覆盖 范围 为 
0.1~400 MeV, 共 20 个 固定 网 格 点 , 近似 于 能 量 对 数 
均匀 网 格 点 . LO 网 格 点 覆盖 范围 为 1.14~6.6, 采用 非 
固定 非 均匀 网 格 点 , 不 同 能 量 其 L 网 格 覆 盖 范 围 不 
同 , 网 格 点 数目 也 不 一 样 , 不 过 在 工大 于 1.2 的 区 域 
P, L 网 格 点 的 位 置 基本 相同 且 间 隔 均 约 为 0.1. B/ 
Bo 网 格 点 覆盖 范围 为 1.0 至 损失 锥 角 对 应 的 B/ Bo 
E, 采用 的 也 是 非 固定 非 均 匀 网 格 点 , 而 且 不 同 工 值 
的 损失 锥 角 不 同 , 其 网 格 点 覆盖 范围 及 网 格 点 设置 也 
不 同 . 同一 根 磁力 线 上 某 一 能 量 的 B/ Bo 网 格 点 , f 
照 以 10 为 底 的 对 数 (全 向 积分 通 量 ) 以 0.25 值 衰减 
设置 , 直到 其 对 应 的 损失 锥 角 所 对 应 的 B/Bo. 详细 
的 AP-8 模式 后 台 网 格 数 据 库 介绍 可 参考 文献 [2]. 

这 里 采用 了 AP-8 模式 数据 库 中 磁 赤 道 面 上 的 
数据 , 根据 其 全 向 积分 通 量 求 得 全 向 微分 通 量 , 用 于 
中 心 区 域 建 模 及 结果 比较 . 为 尽 可 能 减少 插值 带 来 的 
误差 , 磁 赤 道 面 上 所 采用 的 (BEB,，L，B/Bo) 网 格 点 
与 AP-8 模式 保持 一 致 , 取 其 中 B/Bo = 1 的 数据 . 
对 于 某 个 (Eo, Lo) 网 格 点 , 采用 
F(E, Lo) — F(Eo, Lo) 

E — Eg 

求 得 对 应 的 全 向 微分 通 量 J (Ero, Lo), 其 中 F(Eo, Lo) 
为 由 AP-8 模式 数据 库 所 得 全 疝 积分 通 量 ，E 满 
Æ In(E)— In(Eo) = 0.01, 对 于 最 后 一 个 Eo 网 格 点 
满足 I(E)—- In(Eo) = —0.01, F(E, Lo) 采用 对 数 积 
分 通 量 与 能 量 线性 插值 求 得 , 即 对 数 线性 插值 法 , 这 
与 IRBEM 软件 包 中 onera_desp_lib_Win32._x86.dll 处 
理 方式 一 致 *. 需要 注意 的 是 , 如 果 采 用 双 对 数 线性 
插值 方法 求解 F(Eo, Lo), 其 结果 与 上 述 对 数 线 性 插 
值 方 法 所 得 结果 存在 差异 .这 也 是 数据 库 建 模 的 一 
Ar BE. 
1.2 CRRESPRO 模式 数据 库 

CRRESPRO 模式 包含 CRRESPRO Quiet Hf 
平静 子 模式 和 CRRESPRO Active WREE THAR. 
这 两 个 子 模 式 均 为 静态 平均 模式 , 建 模 方法 与 AP-8 
模式 一 样 , 都 采用 后 台数 据 库 形 式 , 所 不 同 的 是 CR- 


J(Eo, Lo) — 
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RESPRO 模式 采用 了 固定 的 E (22 个 网 格 点 ) 工 (90 
个 网 格 点 ) 和 B/ Bo (34 个 格 段 , BE 35 个 网 格 点 ) 网 
格 点 , 数据 库 保存 的 是 这 些 网 格 点 上 的 全 向 微分 通 量 
值 . CRRESPRO 后 台 网 格 数据 库 E. 网 格 点 覆盖 范 
围 为 1.5~81.3 MeV, 近似 能 量 对 数 均匀 网 格 ; L 网 格 
点 覆盖 范围 为 1.0~5.45, 采用 间隔 为 0.05 的 均匀 网 
格 ; B/ Bo 网 格 点 覆盖 范围 为 1.0~684.6, 采用 非 均匀 
网 格 . 详细 的 CRRESPRO 模式 后 台 网 格 数据 库 介 
绍 可 参考 文献 [8]. 
13 RBSP 卫星 REPT 高 能 质子 观测 数据 

辐射 带 风 暴 探 测 器 (Radiation Belt Storm 
Probes, RBSP) 后 改名 为 范 艾 伦 探测 器 (Van Allen 
Probes). RBSP 共 包 含 A, B 两 颗 挛 生 卫 星 , 于 2012 
年 8 H 30 日 发 射 ,采用 一 箭 双 星 方式 送 入 预定 轨 
道 . 这 两 颗 卫 星 轨 道 为 大 椭圆 轨道 , A 卫星 轨道 远 
地 点 为 30050km、 近 地 点 500km; B 卫星 轨道 远 地 
点 为 31250km、 近 地 点 675km, 两 卫星 轨道 倾角 约 
为 10°, 轨道 周期 约 9hP3 

每 颗 RBSP 卫星 上 均 搭载 有 5 种 探测 仪器 , 其 
中 高 能 粒子 、 组 成 与 热 等 离子 体 仪器 (Energetic Par- 
ticle, Composition, and Thermal Plasma, ECT)?! 
和 相对 论 质 谱 仪 (Relativistic Proton Spectrometer, 
RPS) 可 提供 高 能 质子 观测 数据 ，ECT 中 相对 论 
电子 质子 望远镜 (Relativistic Electron Proton Tele- 
scope, REPT)P4 提供 21.25~208.00MeV 高 能 质子 
观测 数据 , RPS 提供 60 MeV~2 GeV 高 能 质子 观测 
数据 , 不 过 后 者 的 数据 目前 尚未 完全 发 布 . 

这 里 将 用 到 RBSP A 卫星 上 2012 年 9 月 1 
日 至 2014 ^E 2 月 22 日 期 间 REPT 的 高 能 质子 
观测 数据 (RBSP B 卫星 观测 结果 类 似 ). REPT 2H 
远 镜 共 提供 8 个 能 道 高 能 质子 方向 微分 通 量 观 测 
数据 ， 分 别 为 21.25, 27.6, 35.9, 46.7, 60.7, 78.9, 
102.6, 208 MeV, HFI AB 77 P8] BADULUU FS TIE E] E 
实现 , 自 旋 周 期 为 10~15sP20. 目前 已 公布 的 最 新 数 
据 为 rbspa.rel02.ect-rept-sci-L3 数据 *, 包括 投掷 角 
归 组 处 理 后 的 数据 和 未 归 组 处 理 的 原始 数据 ， 这 里 
采用 的 是 后 者 中 的 FPDU.Unbinned.0to180 原始 数 
ju. 在 该 数据 中 , 卫星 每 自 旋 一 周 ， 完成 一 个 时 刻 
的 整个 投掷 角 方 向 探测 ， 在 每 一 时 刻 的 数据 中 收集 
了 36 个 投掷 角 方向 数据 , 由 于 卫星 自转 轴 方 向 几乎 


*hppt://www.rbsp-ect.lanl.gov/data.pub/rbspa/rept/level3/ 
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与 磁力 线 垂 直 ， 所 以 每 一 时 刻 基本 可 以 提供 近乎 覆 
TE 0° ~ 180? 的 投掷 角 观 测 数据 , 不 过 不 同时 刻 其 36 
个 投掷 角 的 方向 是 不 一 样 的 . 

在 本 研究 中 , 只 需要 磁 赤 道 面 上 各 能 量 E MOL 
值 高 能 质子 全 向 微分 通 量 的 平均 结果 , 因此 对 原始 数 
据 进行 了 如 下 处 理 . 

(1) 去 除 2012 年 9 月 1 日 至 2014 年 2 月 22 日 
期 间 发 生 太阳 质子 事件 的 观测 数据 . 

(2) 提取 出 B/Bo < 1.001 的 数据 , 作为 磁 赤 道 
面 上 的 观测 结果 . 

(3) 将 磁 赤道 面 上 的 数据 按 工 值 进行 归 组 , 工 履 
盖 范 围 为 1.00~6.20, 分 组 间隔 为 0.01. 

(4) 对 于 各 个 能 道 , 将 每 组 L 中 的 数据 按 投掷 
角 进 行 分 组 , 投掷 角 覆 盖 范 围 为 1? 1797, 分 组 间隔 
为 1°, 然后 将 每 个 投 搓 角 组 中 的 数据 按 观 测 所 得 单 
向 微分 通 量 大 小 进行 排列 , 各 滤 除 上 下 两 端 2596 的 
数据 , 以 确保 尽 可 能 去 除 污 染 数据 , 同时 保留 无 污染 
观测 数据 的 分 布 特性 . 

(5) 对 于 各 个 能 道 , 将 该 工 组 中 所 得 各 投 搓 角 组 
内 滤 除 后 剩 下 的 50% 数据 求 平均 值 , 该 值 即 为 该 L 
值 和 投掷 角 组 中 所 得 平均 单 向 微分 通 量 , 将 用 于 求解 
全 向 微分 通 量 . 

(6) 根据 该 L 组 中 所 得 的 各 投掷 角 平 均 单 向 微 
分 通 量 数据 , 采用 jla) = ji (sino)" +c WE RZ, 
利用 Levenberg-Marquardt 非 线 性 拟 合 方法 , 获得 拟 
合 参 数 jL, n 和 c 的 值 , 其 中 jj 为 垂直 微分 通 量 , a 
为 投掷 角 , (0) 为 a 投掷 角 方向 上 的 微分 通 量 , c 可 
认为 是 数据 平均 污染 值 . 

(7) 利用 


27 T 
J- | | ji (sin a)” sin adad0 = 
0—0 Ja —0 
n/2 
an| jı (sin a)"*!do, 
0 


获得 该 工 值 上 的 全 向 微分 通 量 . 

(8) 其 余 各 个 工 组 也 采用 类 似 方法 处 理 . 

图 1 给 出 了 磁 赤 道 面 上 E = 208 MeV, L= 1.39 
时 投掷 角 数 据 分 布 . 图 1(a) 中 绿色 “+” 为 经 过 以 
上 (1)~(3) 步骤 后 未 经 处 理 的 投掷 角 数 据 分 布 ， 从 
图 1(a) 中 绿色 “+” 分 布 可 以 看 出 , 有 相当 一 部 分 观 
测 数 据 异常 . 图 1(b) 中 绿色 “+” 为 经 过 以 上 (1)~(4) 
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图 1 磁 赤 道 面 上 RBSP A 卫星 观测 所 得 E = 208 MeV, L = 1.39 组 投掷 角 数 据 分 布 . (a) 未 经 处 理 的 投掷 角 数 据 分 布 ， 
(b) 处 理 后 的 投掷 角 数据 分 布 及 拟 合 结果 
Fig.1 Pitch angle distribution of the bin with E = 208 MeV and L = 1.39 from RBSP A observation on the 


magnetic equator: (a) original data, (b) processed data and fitting result 


步骤 处 理 后 的 投掷 角 数据 分 布 , 可 以 看 出 , 基本 滤 除 
了 所 有 有 蜡 常 投掷 角 观 测 数据 , 投 撞 角 数 据 分 布 也 相当 
有 规律 . 图 1(b) 中 黑色 “+” 为 经 过 步骤 (1)~(5) 处 
理 后 所 得 的 各 投掷 角 平 均 单 向 微分 通 量 . 图 1(p) 中 
红线 为 黑色 “+” 点 经 过 步骤 (1)~(6) 后 所 得 的 拟 合 
曲线 结果 . 从 图 中 可 知 , j(Q) = fi (sin a)” +e WEA 
效果 还 是 相当 好 的 . 图 1(b) 中 两 条 蓝 色 竖 线 为 损失 
锥 角 , 该 损失 锥 角 是 在 假定 地 磁场 为 偶 极 场 且 高 能 质 
子 只 能 运行 到 100km 高 度 条 件 下 获得 的 . 从 图 1 可 
知 , 损失 锥 内 实际 观测 的 结果 并 不 为 0, 可 能 是 由 于 
探测 需 受 到 污染 引起 的 , 但 图 中 所 显示 的 损失 锥 内 污 
染 还 比较 小 , 同时 污染 随 投 搓 角 变化 也 较 小 . 通过 步 
又 (7) 既 可 去 除 部 分 污染 , 还 可 求 得 全 向 积分 通 量 . 
以 上 所 用 的 投 撞 角 数 据 、B/Bo 值 、 工 值 均 由 rb- 
spa_rel02_ect-rept-sci-L3 数据 提供 , 通过 OP77Q 外 
源 场 模式 及 IGRF 内 源 场 模式 计算 得 到 95-29. 


在 REPT 观测 所 得 8 个 能 道 高 能 质子 观测 数据 
中 , 根据 以 上 方式 处 理 可 知 , 前 5 个 能 道 损失 锥 内 所 
受 污 染 比 较 严 重 , 而 且 在 相当 一 部 分 区 域 , 污染 随 投 
掷 角 变 化 比较 剧烈 . 因此 本 研究 中 只 采用 78.9, 102.6 
和 208 MeV 三 个 能 道 的 数据 结果 . 


2 iiL AP-8 与 CRRESPRO 
模式 结果 比较 


关于 AP-8 模式 与 CRRESPRO 模式 的 结果 比 
较 以 及 两 模式 所 得 结果 与 卫星 实测 数据 之 间 的 比较 ， 
已 有 相关 文献 研究 02:0. 不 过 其 主要 是 集中 于 模式 
特点 、 南 大 西洋 异常 区 积分 通 量 或 典型 轨道 上 能 谱 
计算 结果 比较 . 这 里 将 集中 于 磁 赤 道 面 上 两 模式 各 能 
段 微分 通 量 观 测 结果 的 差异 比较 . 为 尽 可 能 避免 网 格 
插值 带 来 的 误差 (不 同 的 插值 方法 获得 的 结果 是 有 
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差异 的 ), 保持 模式 原先 的 ERI L 网 格 点 , 只 有 在 利 
用 AP-8 模式 求解 全 向 微分 通 量 时 采用 了 插值 处 理 ， 
如 1.1 节 所 述 . 

图 2(a) 中 显示 了 磁 赤 道 面 上 AP-8 模式 所 得 全 
向 微分 通 量 7 的 空间 分 布 ， 在 图 2(a) F, 黑 线 表 
示 AP-8 MIN 子 模式 各 网 格 点 上 所 得 的 全 向 微分 通 


1019 
10? 


10* 
107 
106 
105 
10* 
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量 结果 , 其 中 网 格 点 如 黑色 “+” 所 示 , 在 绝 大 部 分 区 
域 中 , 其 工 值 的 网 格 点 间隔 为 0.1; 红线 表示 AP-8 
MAX 子 模式 各 网 格 点 上 所 得 的 全 向 微分 通 量 结果 ， 
其 网 格 点 绝 大 部 分 与 AP-8 MIN 子 模 式 一 致 , 网 格 
点 差异 主要 出 现在 工科 146 KEE; 图 中 数字 表示 
各 L-J 曲线 所 对 应 的 能 量 值 . 从 图 2(a) 中 可 知 , 在 


(a) 
AP-8 MIN & MAX Models 


a 
o 
= 
Ins 
d 
E 
E 
E 
t CRRESPRO Quiet & Active Models 
E E- (1.5, 2.1, 2.5, 2.9, 3.6, 4.3, 5.7, 6.8, 8.5, 9.7, 10.7, 132, 
E 16.9, 19.4, 26.3, 30.9, 36.3, 41.1, 47.0, 55.0, 65.7, 81.3} 
图 2 min L AP-8 模式 与 CRRESPRO 模式 所 得 全 向 微分 通 量 空间 分 布 : (a) AP-8 MIN (H2) 和 AP-8 MAX 
(红线 ) 子 模式 所 得 结果 , (b) CRRESPRO Quiet ( 黑 线 ) 和 CRRESPRO Active (红线 和 蓝 线 ) 子 模式 所 得 结果 


Fig.2 Spatial distributions of omni-directional proton differential flux on the magnetic equator from AP-8 and 
CRRESPRO models: (a) AP-8 MIN (black curves) and AP-8 MAX (red curves) submodels; 
(b) CRRESPRO Quiet (black curves) and CRRESPRO Active (red and blue curves) submodels 
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绝 大 部 分 磁 赤 道 面 空间 区 域 ，AP-8 MIN FHA 
和 AP-8 MAX 子 模式 所 得 的 全 向 微分 通 量 结果 基本 
相同 , 其 差异 主要 出 现在 L < 1.3 的 区 域 范围 内 ; 各 
个 能 道 辐射 通 量 的 工 空间 分 布 范 围 基本 上 随 着 能 量 
的 增加 而 减 小 , 每 个 能 道 辐射 中 心 区域 L 值 随 着 能 
量 增加 而 减少 , 辐射 全 向 微分 通 量 也 随 着 能 量 增 加 而 
减少 . 

图 2(b) 中 显示 了 磁 赤 道 面 上 CRRESPRO 模式 
所 得 全 向 微分 通 量 的 空间 分 布 . 在 图 2(b) F, 黑 线 表 
示 CRRESPRO Quiet 子 模式 各 网 格 点 上 所 得 的 全 向 
微分 通 量 结果 , 红线 表示 CRRESPRO Active 子 模式 
各 网 格 点 上 的 全 向 微分 通 量 结果 , 蓝 线 表示 CRRE- 
SPRO Active 子 模式 在 E — 26.3 MeV 时 的 全 向 微分 
通 量 , 两 个 子 模式 的 网 格 点 工 间隔 均 为 0.05; 图 中 各 
曲线 对 应 能 量 见 图 注 , 全 向 微分 通 量 最 强 的 曲线 
对 应 能 量 最 小 , 随 着 能 量 的 增加 , 其 曲线 的 全 向 微分 
通 量 基本 是 不 断 减 小 的 (在 第 二 质子 辐射 带 有 些 异 
9$). 从 图 2(b) 可 知 , 在 磁 赤 道 面 上 绝 大 部 分 的 区 域 
空间 , CRRESPRO Quiet 子 模式 与 CRRESPRO Ac- 
tive 子 模式 所 得 的 全 向 微分 通 量 结果 差异 还 是 比较 
大 的 ; 在 E 2 26.3 MeV 能 道上 , CRRESPRO Active 
子 模式 得 到 的 全 向 微分 通 量 随 LL 变化 的 曲线 出 现 明 
显 的 双 峰 结构 , 即 出 现 了 第 二 质子 辐射 带 , 第 二 峰 所 
在 位 置 大 致 在 L 为 2.22.3 pj, 然而 CRRESPRO 
Quiet. 子 模式 没有 得 到 双 峰 结构 , 但 在 L > 1.9 区 
域 上 , CRRESPRO Quiet 子 模式 得 到 全 向 微分 通 量 
随 着 L 衰减 的 速度 明显 减缓 , 其 结构 似乎 可 以 看 成 
是 原 质子 辐射 带 与 隐 性 第 二 质子 辐射 带 徐 加 造成 的 ， 
只 是 隐 性 的 第 二 质子 辐射 带 峰 值 通 量 明显 比 原 质子 
辐射 带 的 峰值 通 量 小 ,因此 未 能 显示 出 显 性 的 双 峰 
结构 ; 在 E < 26.3MeV 能 道上 , CRRESPRO Qui- 
et 子 模式 和 CRRESPRO Active 子 模式 同样 得 到 隐 
性 的 第 二 质子 辐射 带 结 构 ， 有 些 能 道上 其 至 已 得 到 
显 性 的 双 峰 结构 , 例如 CRRESPRO Quiet 子 模式 所 
得 E — 6.8MeV 能 道上 全 向 微分 通 量 两 个 峰值 大 
小 接近 , 只 是 相距 较 近 , CRRESPRO Active FR 
Æ E = 10.7 MeV 能 道上 也 有 类 似 结 果 ; 对 于 CRRE- 
SPRO Quiet 子 模式 , TE E < 6.8 MeV 能 道上 , 隐 性 的 
质子 辐射 带 出 现在 内 侧 , 在 E > 6.8MeV 能 道上 , ES 
性 的 质子 辐射 带 出 现在 外 侧 ; 对 于 CRRESPRO Ac- 
tive 子 模式 来 讲 , 在 E < 10.7 MeV 能 道上 , 隐 性 的 
质子 辐射 带 出 现在 内 侧 , 在 10.7 < E < 26.3 MeV 能 
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道上 , 隐 性 的 质子 辐射 带 出 现在 外 侧 ; CRRESPRO 
Quiet 与 CRRESPRO Active 子 模式 最 大 的 差异 出 
现在 外 侧 的 质子 辐射 带 上 ; TE L « 1.65 KE, 各 能 道 
上 两 子 模式 所 得 全 向 微分 通 量 差异 较 小 , 对 于 能 量 较 
大 的 能 道 (如 E > 26.3MeV), Æ 1.65 < L < 1.9 IX 
E, 两 子 模式 所 得 全 向 微分 通 量 差异 也 是 比较 小 的 . 
从 图 2(a) 与 图 2(b) 的 比较 可 知 , 从 基本 结构 上 
看 , AP-8 模式 和 CRRESPRO 模式 所 得 全 向 微分 通 
量 差 异 还 是 比较 明显 的 , 主要 表现 为 , 对 于 上 > 2 的 
所 有 能 道 , AP-8 模式 只 表现 出 一 个 质子 辐射 带 的 结 
构 特 征 , 然而 对 于 CRRESPRO Quiet 模式 , 当 能 量 
大 于 一 定 值 时 , 似乎 表现 出 两 个 质子 辐射 带 的 苔 加 结 
构 , 对 于 CRRESPRO Active 模式 , 当 E > 26.3 MeV 
时 , 出 现 了 明显 的 双 峰 结构 ; 在 工大 于 峰值 的 区 域 范 
围 内 , 对 于 各 能 道 , AP-8 模式 得 到 的 对 数 全 向 微分 
MEE L 值 的 增加 先 减 小 得 较 慢 , 而 后 减 小 得 较 快 ， 
即 斜率 越 来 越 大, 然而 CRRESPRO 模式 得 到 对 数 全 
向 微分 通 量 过 了 峰值 一 小 段 区 域 后 随 着 工 值 的 增加 
基本 呈 线 性 方式 减 小 . 如 果 要 定量 比较 , 工 值 较 小 时 ， 
CRRESPRO 模式 所 得 结果 明显 比 AP-8 模式 所 得 结 
ARAS 在 工 值 较 大 时 , CRRESPRO 模式 所 得 结果 明 
显 比 AP-8 模式 所 得 结果 小 (第 二 质子 辐射 带 峰 值 及 
其 附近 区 域 除外 )， 对 于 定量 比较 , 不 同 能 道上 的 比 
较 结果 有 各 目的 差异 表现 , 在 此 不 一 一 讨论 . 第 3 8 
将 重点 比较 AP-8 模式 与 CRRESPRO 模式 所 得 各 
能 道上 全 向 微分 通 量 辐射 中 心 区 域 定量 结果 的 差异 . 


3 质子 辐射 带 中 心 区 域 建 模 


从 第 2 节 可 知 ，AP-8 模式 与 CRRESPRO 模 
式 所 得 全 向 微分 通 量 结果 的 差异 比较 大 .这 里 将 利 
用 AP-8 模式 数据 库 建立 函数 形式 的 质子 辐射 带 中 
心 区 域 模式 (最 强 区 域 在 B/ Bo — 1.0 的 磁 赤 道 面 上 )， 
包括 各 能 道 全 向 微分 通 量 峰值 所 对 应 的 工 值 (Le) 随 
能 量 己 的 变化 模式 及 各 能 道 全 向 微分 通 量 峰 值 Jmax 
随 能 量 E 的 变化 模式 . 在 构建 辐射 中 心 区 域 模式 过 
程 中 , 定量 比较 AP-8 模式 与 CRRESPRO 模式 所 得 
质子 辐射 带 辐射 中 心 区 域 结 果 的 差异 . 
3.1 ”辐射 中 心 区 域 L. 随 E 的 变化 模式 

根据 图 2 所 得 结果 , 容易 获取 AP-8 模式 和 CR- 
RESPRO 模式 所 得 磁 赤 道 面 上 各 能 道 高 能 质子 最 强 
全 向 微分 通 量 所 对 应 的 LB (Lc) 及 其 相应 全 向 微 
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分 通 量 峰值 .sx. 对 于 CRRESPRO Active 模式 , 在 
高 能 段 出 现 两 个 峰值 , 只 采用 内 侧 较 大 的 峰值 及 其 所 
对 应 的 Le fH. 

图 3 显示 了 AP-8 模式 和 CRRESPRO 模式 所 得 
磁 赤 道 面 上 各 能 道 高 能 质子 最 强 全 向 微分 通 量 所 对 
应 的 工 值 随 能 量 五 的 变化 . 从 图 中 AP-8 模式 (AP-8 
MIN 与 AP-8 MAX 子 模式 所 得 结果 一 致 ) 所 得 结果 
T4, Æ E < 1MeV Ñ E >20MeV Kit, 辐射 带 中 
心 位 置 Lo 随 能 量变 化 相对 比较 缓慢 ; 在 低能 端 区 域 ， 
L. 随 着 能 量 减 小 其 增加 的 速度 越 来 越 缓慢 , 有 趋 近 
于 饱和 的 趋势 ; 在 高 能 端 区 域 , Lo 随 能 量 的 增加 而 减 
小 的 速度 越 来 越 缓慢 ,也 有 趋 近 于 饱和 状态 的 趋势 ; 
Æ 1 < E « 20MeV Kik, 辐射 带 辐射 中 心 位 置 随 
着 能 量 增 加 而 减 小 的 速度 很 快 , L. 与 lg(B) 基本 上 
呈 线 性 变化 关系 . 从 图 中 CRRESPRO 模式 所 得 结 
果 可 知 , Quiet FRIE E X 5.7~6.8 MeV EF Le H 
现 明 显 突变 ; Active 子 模式 在 E 为 10.7—13.2MeV 
时 Ze 也 出 现 突变 情况 ; 在 突变 所 对 应 的 能 量 之 前 ， 
CRRESPRO 模式 也 显示 出 Le 与 lg(B) 基本 呈 线 性 
变化 关系 , 在 突变 对 应 的 能 量 之 后 , Ze 随 着 能 量 增加 
而 减 小 的 速度 相对 比较 缓慢 . CRRESPRO 模式 所 得 
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AP-8 MIN: © 

AP-8 MAX: 十 
CRRESPRO Quiet: ¥ 
CRRESPRO Active: 4- 
Fit: 
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到 的 这 种 突变 现象 主要 是 由 于 隐 性 的 第 二 质子 辐射 
带 和 原 质子 辐射 带 炙 加 造成 的 , 如 图 2(b) 所 示 . 

由 于 AP-8 模式 能 量 的 覆盖 范围 较 广 , 而 且 Le 
与 五 之 间 比 较 有 规律 性 , 这 里 采用 AP-8 模式 数据 
库 所 得 结果 进行 辐射 带 中 心 区 域 Le 5j E BEES. 根 
据 图 3 中 AP-8 模式 所 得 Le 和 E 的 关系 , 选用 

a 
LE) = DIR 十 as 

MERE, KH Levenberg-Marquardt 非 线性 拟 合 
方法 , 得 到 各 拟 合 系数 分 别 为 ao = 一 0.515, al = 
—1.195, a3 = 0.245, a3 = 3.366, Le 拟 合 的 标准 偏 
2824 0.0850. 拟 合 函 数 所 得 到 的 曲线 如 图 3 中 黑色 
曲线 所 示 , AP-8 模式 结果 与 拟 合 结果 吻合 很 好 , AP- 
8 模式 结果 基本 均匀 分 布 于 拟 合 函 数 所 得 曲线 的 两 
边 , 而 且 Lc 拟 合 的 标准 偏差 比 AP-8 模式 辐射 中 心 
及 其 附近 L 网 格 点 间隔 0.1 还 要 小 . 所 得 Le 和 五 
的 关系 模式 修复 了 AP-8 模式 因 网 格 点 粗糙 所 带 来 
的 小 幅度 突变 , 同时 模式 也 具备 一 定 的 适用 范围 
拓展 能 力 . 通过 CRRESPRO 模式 结果 与 拟 合 函数 
曲线 的 比较 可 知 , 在 大 部 分 能 道上 CRRESPRO 模 
式 所 得 Le 比 函数 拟 合 模式 所 得 结果 小 . 


10! 10? 10? 10* 
E/MeV 


图 3 AP-8 模式 和 CRRESPRO 模式 所 得 磁 赤 道 面 上 各 能 道 高 能 质子 最 强 全 向 微分 通 量 对 应 工 值 (Lc) 
随 能 量 E 的 变化 以 及 Lc(E) 函数 拟 合 结果 


Fig.3 Distributions of Le corresponding to the peak value of omni-directional proton differential flux 


for various energy channels in the E-L coordinates on the magnetic equator from AP-8 and 
CRRESPRO models and the fitting result of L.(E) function 


HEF: 质子 辐射 带 辐射 中 心 区 域 模型 


3.2 ”辐射 中 心 区 域 Jma 随 E 的 变化 模式 

图 4 显示 了 AP-8 模式 和 CRRESPRO 模式 所 得 
磁 示 道 面 上 各 能 道 高 能 质子 最 强 全 向 微分 通 量 Jmax 
随 能 量 E 的 变化 .从 图 中 AP-8 模式 (AP-8 MIN 
和 AP-8 MAX 子 模式 所 得 结果 一 致 ) 的 结果 可 以 
看 出 ， 随 着 能 量 增加 ， 辐 射 带 中 心 区 域 最 强 全 向 微 
分 通 量 Jma 下 降 速 度 非常 快 , 4 E = 0.01 MeV 
Hf, Jmax = 1.43 x 10?.cm-7?7.s-|. MeV, 4 E = 
400MeV 时 , Jmax = 6.94cm7?:s-1-MeV-!; Æ E < 
50 MeV 区 域 , 1gJ 与 lgE 呈 类 似 抛 物 线 关系 ,在 已 之 
50 MeV 区 域 , 二 者 呈现 出 另 一 个 类 似 抛 物 线 关系 . 从 
图 中 CRRESPRO 模式 所 得 结果 可 知 , lg(J) 与 lg(E) 
呈现 出 近似 线性 关系 . 通过 AP-8 模式 与 CRRE- 
SPRO 模式 所 得 结果 比较 可 知 , 在 E « 10 MeV 区 间 
内 , CRRESPRO Quiet 子 模式 所 得 Jmax 明显 比 AP- 
8 模式 所 得 结果 小 ，CRRESPRO Active 子 模式 所 
得 Jmax 大 部 分 与 AP-8 模式 所 得 结果 一 致 ; 在 E > 
10 MeV 区 间 , CRRESPRO Quiet 子 模式 与 CRRE- 
SPRO Active 子 模式 所 得 Jaax 明显 比 AP-8 模式 所 
得 结果 大 . 
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仍然 采用 AP-8 模式 数据 库 所 得 结果 进行 辐射 
带 辐 射 中心 区 域 Jma 与 E 的 关系 建 模 . 根据 图 4 
中 AP-8 模式 所 得 Jua. 与 E 的 关系 , 选用 


In(Jmax) = bo + bie? F" +b4E 


拟 合 方程 , 采用 Levenberg-Marquardt 非 线性 拟 合 方 
ik. 得 到 拟 合 系数 分 别 为 bo = 5.496, bı = 16.858, 
bz = —0.284, bs = 0.592, b4 = 一 0.00865, In(Jmax) 拟 
合 的 标准 偏差 为 0.226. 拟 合 函数 结果 如 图 4 黑色 曲 
线 所 示 , AP-8 模式 结果 与 拟 合 结果 吻合 很 好 , 同时 模 
式 也 具备 一 定 的 值 适用 范围 拓展 能 力 . 


4 ” 磁 赤 道 面 上 RBSP A 卫星 观测 
结果 与 模式 计算 结果 比较 


根据 1.3 节 所 描述 的 处 理 方法 , 容易 获得 RBSP 
A 卫星 上 REPT 望远镜 在 磁 赤 道 面 的 78.9, 102.6 
和 208 MeV 三 个 能 道观 测 所 得 全 向 微分 通 量 随 着 L 
值 的 变化 , 其 结果 如 图 5 中 “十 ”所 示 . 图 5 中 同时 给 
出 了 AP-8 模式 和 CRRESPRO 模式 在 这 三 个 能 道 
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图 4 AP-8 模式 和 CRRESPRO 模式 所 得 磁 赤 道 面 上 各 能 道 高 能 质子 最 强 微 分 通 量 Jmax 
随 能 量 E 的 变化 及 Jaax (I7). 函数 拟 合 结果 


Fig.4 Distributions of the peak value of omni-directional proton differential flux Jmax 


for various energy channels in the E-J coordinates on the magnetic equator from AP-8 and 


CRRESPRO models and the fitting result of Jmax(E) function 


302 


E=78.9 MeV 


-2.g-l -MeV-!) 


differential flux/{ 


E=208 MeV 


1.0 135 2.0 
E 


E=102.6 MeV 


Chin. J. Space Sci. 空间 科学 学 报 2015, 35(3) 


(a) 
RBSP A: + 
CRRESPRO Active: 
CRRESPRO Active: ——— 
AP-8 MIN: ——— 
AP-8 MAX: ——— 


RBSP A: + 
AP-8 MIN: 


RBSP A: + 
AP-8 MIN: 
AP-8 MAX: ——— 


n 
* 
为 
2 
+4 


十 
+ 


2.5 3.0 3:5 


图 5 磁 赤 道 面 上 RBSP A 卫星 REPT 望远镜 三 个 能 道观 测 所 得 全 向 微分 通 量 结果 与 AP-8 模式 和 


CRRESPRO 模式 计算 所 得 结果 比较 


Fig.5 Comparison of omni-directional proton differential flux from the observation of RBSP A satellite's REPT 


telescope and from the calculations of AP-8 and CRRESPRO models for the three energy channels 


上 的 计算 结果 由 于 这 三 个 能 道 并 不 在 AP-8 模式 
和 CRRESPRO 模式 数据 能 基 网 格 点 上 , 因此 采用 
对 数 全 向 微分 通 量 与 对 数 能 量 的 线性 插值 求 得 ， 即 
双 对 数 插值 方法 .。 对 于 AP-8 MIN 模式 或 AP-8 
MAX 模式 , 不 同 能 道 其 L 网 格 点 在 L 值 较 小 时 是 
不 同 的 , 因此 先 采 用 对 数 全 向 微分 通 量 与 工 值 的 线 
性 插值 求 得 , 获得 上 下 两 能 道 同样 L 网 格 点 数据 , 再 
采用 双 对 数 线 性 插值 方法 求 得 磁 赤 道 面 78.9, 102.6 


或 208 MeV 能 道上 模式 计算 所 得 全 向 微分 通 量 随 L 
值 变 化 结果 . 由 于 CRRESPRO 模式 能 量 使 用 范围 
ERY 81.3 MeV, 仅 能 获得 78.9 MeV 能 道 结果 , 如 
图 5 所 示 . 

从 图 5 可 以 看 出 , 在 磁 赤 道 面 上 RBSP A E 
星 观 测 所 得 78.9, 102.6 和 208 MeV 三 个 能 道上 的 
观测 结果 明显 比 AP-8 MIN 模式 或 AP-8 MAX 模 
式 计 算 结果 偏 大 约 几 倍 到 十 几 倍 ; E L < 1.2 KI, 


陶 宏 任 等 : 质子 辐射 带 辐 射 中 心 区 域 模型 


RBSP A 卫星 观测 所 得 78.9 MeV 能 道上 的 观测 结果 
与 CRRESPRO Quiet 子 模式 或 CRRESPRO Active 
子 模 式 结果 相当 , 而 其 他 区 域 也 明显 比 CRRESPRO 
Quiet 子 模式 或 CRRESPRO Active 子 模式 计算 所 
得 结果 大 ; RBSP A 卫星 在 78.9, 102.6 和 208 MeV 
三 个 能 道上 均 观 测 到 隐 性 的 第 二 质子 辐射 带 ， 该 结 
R- CRRESPRO Quiet 模式 所 得 结果 一 致 , 而 且 所 
得 全 向 微分 通 量 随 工 值 衰减 速度 的 明显 减缓 均 出 现 
在 1.9 < L< 2.0 区 域内 . 

从 图 3 可 以 明显 看 出 ，RBSP A 卫星 观测 所 
得 78.9, 102.6 和 208MeV 三 个 能 道上 最 强 全 向 微 
分 通 量 对 应 的 工 值 与 AP-8 模式 所 得 结果 基本 一 致 ， 
比 本 文 所 得 辐射 中 心 区 域 模式 所 得 Ze 略 小 , 差 值 
在 0.05 以 内 , 在 AP-8 模式 此 处 工 网 格 点 间隔 的 0.1 
以 内 . 从 图 4 可 以 明显 看 出 ，RBSP A 卫星 观测 所 
得 78.9, 102.6 和 208 MeV 三 个 能 道上 最 强 全 向 微 
分 通 量 明 显 大 于 AP-8 模式 结果 或 本 文 辐射 中 心 区 
域 模式 结果 , 分 别 为 本 文 辐射 中 心 区 域 模式 计算 结果 
的 5.7 倍 、4.1 倍 和 12.6 fr. 


5 讨论 及 分 析 


利用 AP-8 模式 与 CRRESPRO 模式 数据 库 , 对 
磁 赤 道 面 上 各 能 道 高 能 质子 全 向 微分 通 量 结果 进 行 
了 比较 . AP-8 MIN 和 AP-8 MAX 两 子 模 式 在 绝 
大 部 分 磁 赤 道 面 空间 所 得 全 向 微分 通 量 结 果 基 本 相 
H, 其 差异 主要 出 现在 L < 1.3 的 区 域内 ,如 图 5 
所 示 . 其 原因 为 AP-8 模式 主要 考虑 的 是 太阳 活动 
周 对 低 高 度 质子 辐射 带 的 影响 .在 低 高 度 上 AP-8 
MIN 子 模式 所 得 结果 之 所 以 比 AP-8 MAX 子 模式 
高 ,依据 Miyoshi 等 研究 l, 其 原因 有 两 个 方面 . 
(1) 太阳 活动 低 年 时 产生 质子 辐射 带 的 宇宙 线 中 子 反 
照 增强 , 即 在 太阳 活动 低 年 时 , 太阳 风 和 行星 际 磁场 
较 弱 且 比 较 稳定 , 对 宇宙 线 调 制作 用 较 弱 , 使 得 到 达 
地 球 的 宇宙 线 较 强 , 产生 了 强 的 宇宙 线 中 子 反 照 , 造 
成 低 高 度 上 AP-8 MIN 模式 结果 较 大 ; (2) 太阳 活 
动 高 年 时 质子 辐射 带 损失 增强 ， 即 在 太阳 活动 高 年 
时 , 太阳 EUV 辐射 增强 , 使 得 大 气温 度 升 高 , 大 气 的 
标高 升 高 , 大 气 密度 上 升 , 造成 质子 辐射 带 中 有 更 多 
的 高 能 质子 被 耗 散 到 损失 锥 内 . CRRESPRO Quiet 
和 CRRESPRO Active 两 子 模式 在 绝 大 部 分 磁 赤 道 
面 空间 所 得 的 全 向 微分 通 量 结果 差异 比较 大 , 特别 是 
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1E E 2 26.3 MeV 能 道上 , 从 CRRESPRO Quiet 子 模 
式 所 得 的 隐 性 第 二 质子 辐射 带 发 展 到 CRRESPRO 
Active 子 模式 所 得 的 显 性 第 二 质子 辐射 带 ， 对 于 第 
二 质子 辐射 带 的 形成 , 认为 需要 强 太阳 质子 事件 、 大 
地 磁 暴 和 大 行星 激 波 三 个 条 件 , 目前 国际 上 主要 有 以 
下 两 种 解释 27-29. 第 一 , 太阳 质子 通过 磁 暴 急 始 被 
脉冲 电场 加 速 并 注入 内 磁 层 , 从 而 形成 第 二 质子 辐射 
带 ; 第 二 , 磁 暴 主 相 期 间 渗 透 到 内 磁 层 的 太阳 质子 在 
磁 暴 恢复 相 期 间 被 捕获 , 从 而 形成 第 二 质子 辐射 带 . 

关于 AP-8 模 式 与 CRRESPRO 模式 之 间 的 明显 
差异 , 其 影响 因素 较 多 , 主要 可 能 有 以 下 四 个 方面 原 
因 . (1) 质子 辐射 带 是 动态 的 , 受 多 种 因素 影响 , 例如 
地 磁 暴 、 太 阳 EUV 辐射 等 , 而 两 模式 都 仅 考虑 了 一 
个 因素 , 而 且 所 考虑 的 影响 因素 不 同 , 此 外 只 考虑 了 
静态 的 平均 结果 ; (2) 两 模式 建 模 数 据 不 一 样 , 其 定 标 
结果 可 能 存在 差异 , 而 且 数 据 覆 盖 的 年 代 也 不 一 样 并 
相距 二 三 十 年 , 质子 辐射 带 在 这 样 长 时 间 过 程 中 可 能 
出 现 演变 ; (3) 不 同 卫星 高 能 质子 探测 器 受 污 染 情 况 
不 一 样 , 即使 是 同一 卫星 高 能 质子 探测 器 , 在 不 同 能 
道 、 不 同位 置 或 不 同 投掷 角 方向 受 污 染 的 情况 也 不 
相同 , 因此 难以 确保 所 采用 的 数据 污染 被 去 除 干 净 ; 
(4) 在 质子 辐射 带 建 模 时 , 除 定 标 和 去 污染 外 ,处 理 
方式 不 当 也 会 带 来 建 模 数据 误差 , 例如 CRRES 卫星 
在 近地点 运行 速度 太 快 而 L 网 格 点 距离 (0.05) JA 
所 带 来 的 误差 、L 和 B/Bo 网 格 化 求 平均 带 来 的 误 
差 以 及 磁场 模式 不 同 造成 的 差异 等 . 

鉴于 目前 AP-8 模式 与 CRRESPRO 模式 整体 
结果 差异 较 大 有 旦 网 格 化 建 模 存在 问题 , 采用 AP-8 模 
式 数据 库 数据 建立 了 函数 形式 的 质子 辐射 带 辐 射 中 
心 模式 , 避免 不 同 能 量 之 间 最 强 全 向 微分 通 量 及 其 对 
应 的 工 值 (Ze) 出 现 明显 突变 等 非 物理 现象 , 也 避免 
了 网 格 化 模式 在 求解 非 网 格 点 上 的 结果 时 不 同 插 值 
方法 带 来 的 结果 差异 , 同时 函数 形式 的 模式 还 可 拓展 
能 量 的 适用 范围 . 对 于 图 3 所 示 AP-8 模式 所 得 不 同 
能 量 Le 的 不 连续 现象 , 似乎 是 由 于 建 模 时 网 格 化 过 
于 粗糙 造成 的 , 因为 AP-8 模式 所 得 Ze 与 本 文 辐射 
中 心 模式 所 得 结果 的 差异 基本 在 AP-8 模式 此 处 网 
格 点 间隔 的 0.1 以 内 . 图 3 所 示 CRRESPRO 模式 
所 得 Le 大 多 比 AP-8 模式 和 本 文 辐射 中 心 模式 所 得 
结果 小 , 这 可 能 与 地 磁场 长 时 间 豪 减 引 起 党 层 下 沉 有 
X. CRRESPRO 模式 所 得 Le 随 能 量变 化 出 现 明显 
突变 , 这 是 由 于 CRRESPRO 模式 隐 性 质子 辐射 带 随 
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着 能 量 上 升 从 内 侧 变化 到 外 侧 造成 的 , 如 图 2 所 示 . 
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据 库 由 美国 空军 研究 实验 室 (AFRL) 提供 , 洛斯 阿拉 莫 斯 国 


对 于 中 心 区 域 各 能 其 最 强 的 全 向 微分 通 量 , 本 文 辐射 
中 心 模式 、AP-8 模式 及 CRRESPRO 模式 所 得 结果 
差异 较 小 ， 本 文 辐射 中 心 模式 关于 能 其 合理 的 拓展 
范围 还 需要 其 他 卫星 提供 低能 段 和 高 能 段 数据 进行 
验证 , 目前 本 文 辐射 中 心 模式 所 得 Jmax 在 高 能 段 , 例 
W E 1GeV Kik, 其 结果 可 能 衰减 过 快 . 在 获取 高 
能 段 实测 数据 之 后 (例如 RBSP RPS 观测 数据 完全 
公开 后 ), 只 需 将 模式 拟 合 方程 由 DAE 改 为 BE” 再 
进行 拟 合 . 

通过 各 模式 所 得 与 RBSP A 卫星 REPT 望远镜 
观测 结果 的 比较 可 知 , 至 少 在 78.9, 102.6 或 208 MeV 
能 道上 RBSP A 观测 所 得 全 向 微分 通 量 过 大 , 需要 
重新 定 标 ， 不 过 这 三 个 能 道 最 强 全 向 微分 通 量 所 对 
应 的 工 值 与 AP-8 模式 或 本 文 辐射 中 心 模式 计算 结 
果 差 异 很 小 ， 因 为 对 于 某 一 能 道 高 能 质子 而 言 ， 如 
果 受 污染 或 定 标 不 准 所 得 的 全 向 微分 通 量 观测 结果 
与 实际 结果 呈 线 性 关系 , 只 要 其 缩放 系数 和 常数 因子 
与 工 值 无 关 , 则 观测 所 得 的 最 强 全 向 微分 通 量 对 应 
的 工 值 就 不 会 因为 污染 或 定 标 不 准 而 改变 . 


6 结语 

依据 AP-8 和 CRRESPRO 模式 数据 库 ， 系统 
比较 了 其 在 磁 赤 道 面 上 计算 结果 的 差异 并 给 出 相应 
解释 . 同时 利用 AP-8 模式 数据 库 数据 建立 起 质子 
辐射 带 函 数 形式 的 中 心 区 域 模式 ， 包括 各 能 道 全 向 
微分 通 量 峰值 对 应 L 值 随 能 量 E 的 变化 模式 以 及 
各 能 道 全 向 微分 通 量 峰 值 .hx 随 E WEER. 
利用 RBSP A 卫星 2012 Æ 9 H 1 HÆ 2014 4E 2 
月 22 日 期 间 REPT 望远镜 磁 赤 道 面 上 78.9, 102.6 
和 208 MeV 能 道观 测 所 得 全 向 微分 通 量 数据 , 与 AP- 
8 模式 、CRRESPRO 模式 以 及 本 文 所 得 辐射 中 心 区 
域 模式 计算 结果 进行 了 比较 , 得 到 这 3 个 能 道上 RB- 
SP A 卫星 观测 所 得 各 能 道 微分 通 量 明显 偏 大 , 最 强 
微分 通 量 对 应 的 工 值 与 AP-8 模式 或 本 文 辐射 中 心 
模式 所 得 结果 基本 一 致 . 此 外 , RBSP A 卫星 也 观测 
到 CRRESPRO Quiet 模式 所 得 到 的 隐 和 性 第 二 质子 
辐射 带 结构 . 


致谢 “AP-8 模式 代码 和 模式 数据 库 由 NASA 空间 科学 
数据 中 心 (NSSDC) 提供 , CRRESPRO 模式 代码 和 模式 数 


家 实验 室 (LANL)ECT 科学 控制 中 心 (SOC) 提供 了 数据 下 
载 , RBSP REPT 相关 团队 对 数据 进行 了 预 处 理 . 
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